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separating  the  intermolecular  potential  in  the  proper  way  to  obtain  most  rapid  con¬ 
vergence  of  perturbation  theories  is  discussed.  It  is  concluded  that  the  structure 
of  simple  liquids  is  determined  mostly  by  the  rapidly  varying  parts  of  the  potential. 
Further,  the  ORPA  provides  an  accurate  theory  for  calculating  the  contributions  from 
different  forces  to  the  structure  and  thermodynamics  of  liquids.  / 
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Abstract 


The  optimized  random  phase  approximation  (ORPA)  is  applied  to  the 
calculation  of  the  thermodynamic  properties  and  pair  correlation  function 
of  simple  liquids.  General  formulas  are  presented,  together  with  results 
ror  the  Leonard- Jones  fluid,  which  are  compared  with  Monte  Carlo  and 
molecular  dynamic  results.  Excellent  agreement  Is  obtained  for  the 
entire  single  phase  fluid  region  of  the. phase  diagram,  except  for  the 
critical  region  and  the  very  low  temperature  vapor.  The  theory  con¬ 
verges  especially  rapidly  at  high  densities  or  at  high  temperatures. 

The  importance  of  separating  the  i ntc -molecular  potential  in  the  proper 
way  to  obtain  most  rapid  convergence  of  perturbation  theories  is  dis¬ 
cussed.  It  is  concluded  that  the  structure  of  simple  liquids  is  deter¬ 
mined  mostly  by  the  rapidly  varying  parts  of  the  potential.  Further, 
the  ORPA  provides  an  accurate  theory  for  calculating  the  contributions 
from  different  forces  to  the  structure  and  thermodynamics  of  liquids. 


I.  Introduction 


In  a  recent  discussion  of  the  equilibrium  structure  of  simple 
liquids',  we  i'ltroduced  a  method  for  separating  the  i ntermolecular 
potential  into  two  parts,  an  attractive  part  and  a  repulsive  part,  and 
discussed  the  conceptual  advantages  of  this  separation.  In  particular 
we  concluded  that  the  structure  of  the  Lennard- Jones  liquid  is  deter¬ 
mined  largely  by  the  repulsive  part  of  th3  potential  at  high  density. 

To  test  this  conclusion,  subsequent  Monte  Carlo  calculations2’3 
were  performed  for  hypothetical  fluid  whose  potential  is  just  the  re¬ 
pulsive  part  of  the  Lennard- Jones  potential.  The  differences  between 
the  structure  of  this  fluid  and  the  structure  of  the  Lennard-Jones 
liquid  (at  the  same  temperature  and  density)  clearly  are  due  to  the 
attractive  forces.  At  high  densities,  these  differences  were  indeed 
found  to  oe  small,  but  not  negligible. 

Any  theory  which  neglects  the  effect  of  attractive  forces  on  the 
structure  of  a  liquid  may  be  called  a  "high-temperature  approximation 
(HTA)",  since  it  is  reasonable  to  expect  such  an  approximation  to  be 
valid  at  high  temperatures.  This  type  of  approximation  may  also  be 
applied  in  thermodynamic  perturbation  theory.  For  the  Lennard-Jones 
fluid,  the  HTA  gives  remarkably  accurate  thermodynamic  properties  not 
only  at  high  temperatures  but  also  at  high  densities  for  all  temperatures 
down  to  and  including  the  triple  point  temperature4’2. 

Thus  for  the  Lennard-Jones  liquid,  the  HTA  is  useful  for  all  high 
temperature  strtes  and  for  all  high  density  states.  Over  this  range 
it  provides  a  qualitatively  correct  theory  of  structure  and  a  quanti¬ 
tatively  accurate  theory  of  thermodynamic  properties. 
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To  obtain  a  more  general  theory  of  liquid  structure  and  thermo¬ 
dynamics,  one  which  is  applicable  to  lower  densities  or  temperature  and 
to  other  types  of  fluids,  it  is  necessary  to  have  an  accurate  method 
for  going  beyond  the  HTA,  i.e.  for  calculating  the  effect  of  attrac¬ 
tive  forces  on  the  liquid  structure.  In  this  paper  we  apply  the  re¬ 
cently  developed5  "optimized  random  phase  approximation  (ORPA)"  to 
this  problem.  For  the  Lennard-Jones  liquid,  the  method  is  very  accurate 
for  the  enti re  single-phase  fluid  region  of  the  phase  diagram  (except 
for  the  critical  region  and  the  vapor  far  below  the  critical  tempera¬ 
ture).  The  thermodynamic  results  agree  with  those  of  Monte  Carlo  and 
molecular  dynamics  calculations  to  within  the  estimated  uncertainty 
of  the  latter,  and  the  pair  correlation  functions  and  structure  fac¬ 
tors  also  agree  well  with  those  obtained  from  the  computer  calculations. 
The  theory  Is  very  rapidly  convergent  for  the  Lennard-Jones  liquid  (much 
more  so  than  the  Barker- Henderson  perturbation  theory6),  and  so  we 
expect  that  it  will  be  applicable  to  a  wide  variety  of  liquids,  in¬ 
cluding  those  in  which  the  effect  of  attractions  is  great.  For  exam¬ 
ple,  ORPA  (in  suitably  generalized  form7)  can  be  applied  to  reasonable 
potential  models8  for  liquid  water.  In  addition,  the  ORPA  provides  a 
way  of  understanding  why  the  effect  of  attractions  is  so  small  for 
simple  liquids. 

In  the  next  section,  we  present  some  preliminary  definitions 
and  review  the  basic  Ideas  of  the  HTA  and  the  ORPA.  Section  III  dis¬ 
cusses  the  calculational  procedures  Involved  in  applying  the  theory. 

The  results  for  the  Lennard-Jones  fluid  are  presented  In  Section  IV, 
and  conclusions  are  discussed  In  Section  V. 
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II.  Definitions  and  General  Theory 
A.  Definitions 


The  general  system  under  consideration  here  Is  a  classical  fluid 
of  N  molecules  in  a  volume  V  (the  average  number  density  Is  p  =  N/V) . 
We  shall  assume  that  the  total  potential  energy,  W,  can  be  written  as 
f.  sum  of  spherically  symmetric  pair  potentials,  w: 


W  “  £  W(rlj)  ’ 

*  <J  =  1 

where  r.j  is  the  distance  between  molecules  I  and  J. 

The  Helmholtz  free  energy  for  the  system  is  given  by 


(1) 


-MU  tf=v-  tn  , 


(2) 


where  AA  denotes  the  excess  free  energy  with  respect  to  the  ideal  gas 
at  the  same  temperature,  density  and  volume,  and  fT1  is  Boltzmann's 
constant  times  the  temperature.  The  thermodynamic  properties  are  obtained 
by  differentiating  ft  with  respect  to  p  and  p.  The  correlation  func¬ 
tions  can  be  obtained  by  functional  differentiation.  In  particular, 

g(r)  =  -(2/p’j  [6&6e»<(r)]  =  h(r)+1  ,  (3) 

where  g(  r)  is  the  radial  distribution  function,  and  h(  r)  is  called  the 
two-particle  correlation  function.  The  dimensionless  Fourier  transform 
of  h(  r)  is 

fyk)  =  P  f  dr  h(0  exp(-lk-r)  =  x( k)  -  1  , 

J  ~  ~  ~ 

where  x(k)  is  called  the  structure  factor. 


(4) 


! 
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The  methods  we  describe  in  this  paper  are  applicable  to  classi¬ 
cal  fluids  for  which  the  pair  interaction  contains  a  harshly  repulsive 
part  that  is  short  ranged  (in  addition  to  seme  attractive  parts  and/or 
long  ranged  repulsions)  : 

w(  r)  «  u0(  r)  +  u(  r)  ,  ( 5) 

where  u0( r)  is  the  harsh  repulsion.  We  call  u0( r)  the  reference  poten¬ 
tial  ;  it  is  the  pair  potential  for  the  reference  fluid  in  which  the 
forces  between  the  molecules  are  short  ranged  and  harshly  repulsive. 

The  remainder,  u( r) ,  is  called  the  perturbation  potential. 

In  saying  that  the  potential  Uo  is  "short  ranged  and  harshly  repul¬ 
sive"  we  mean  that  the  Boltzmann  factor  for  u0,  exp[ -j3u0(  r) ]  ,  closely 
resembles  the  Boltzmann  factor  for  a  hard  sphere  interaction  (which  Is 
a  step  function  rising  from  0  to  1  at  a  distance  corresponding  to  the 
hard  sphere  diameter  d) .  In  particular,  exp[-3u0(r)]  must  be  zero  for 
small  r,  unity  for  large  r,  and  it  must  rise  from  essential ly  zero  to 

essentially  unity  over  a  range  of  r  values  which  is  small  compared 

] 

with  the  r  at  which  It  has  the  value  of -jj . 

Whenw(r)  contains  short  ranged,  strongly  repulsive  forces  for 
0<r<r0  only,  we  want  to  choose  Uo(r)  to  have  the  same  forces,  and  so 
u0  and  u  can  be  defined  by  Eq.  (5)  and  the  following: 

du<>(r)/dr  =  dw(r)/dr,  r<ro 

Uo(r)  =  0  r>  r0  .  (6) 

For  this  separation  it  follows  that 

u(  r)  =  constant  =  u(  r0) ,  r  <  r<>  . 

Also  u(r)  Is  continuous  and  more  slowly  varying  than  If  we  had  not 
included  all  the  short  ranged  repulsive  forces  In  uo(r).  An  example 

I 

-  -----  ■  —  - S - : - -  .  - - - -mmm#  | 
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of  this  kind  of  potential  separation  is  depicted  in  Fig.  1  for  the  case 
in  which  w(  r)  is  the  Lennard- Jones  potential  and  ro  is  chosen  te  be  the 
minimum  of  the  potential.  For  this  separation  of  the  Lennard-Jones 
potential,  the  perturbation  contains  all  the  attractive  forces  and  no 
repulsive  forces  and  is  more  slowly  varying  than  Uy.  In  the  most  gen¬ 
eral  situation  u(  r)  may  also  contain  some  slowly  varying  and/or  long 
ranged  repulsive  interactions. 

There  are  two  major  reasons  for  separating  w( r)  into  a  harshly 
repulsive  reference  potential  plus  a  slowly  varying  perturbation.  First, 
a  fluid  in  which  the  interactions  are  short  ranged  and  harshly  repul¬ 
sive  Is  obviously  similar  to  a  hard  sphere  fluid.9  It  is  thus  possible 
to  exploit  this  similarity  and  define  a  new  fluid,  called  the  trial 
system,  in  which  the  pair  interaction  ?s  given  by 


wT(  0  =  urf(  r)  +  uT(  r) 

where 

Uj.(r)  =  u(r),  r >d  , 

u,(r)  =  ”.  r <d 
=  0,  r >d 


(7a) 

(7b) 

(8) 


and  d<ro-  We  call  this  the  trial  system  because  it  is  necessary  to 
"try"  (see  below)  various  different  values  of  d  to  find  the  one  which 
Is  most  closely  related  T.o  the  fluid  with  the  pair  potential  w(  r) . 

For  distances  smaller  than  d,  the  perturbation  in  the  trial  system, 
u^.(  r) ,  can  take  on  any  finite  value.  The  trial  system  differs  from 
the  fluid  of  interest  only  in  that  u<>  has  been  replaced  by  the  hard 
sphere  potential  u^.  This  difference  makes  the  trial  system  more 
amenable  to  statistical  mechanical  analysis  than  the  actual  fluid  of 
Interest.  In  particular,  the  short  ranged  repulsions  of  the  trial 
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system  are  hard  sphere  repulsions,  and  a  great  deal  of  Information 
about  the  hard  sphere  fluid  Is  available  from  the  results  of  computer 
experiments.  Further,  our  technique  for  describing  the  effects  of 
the  perturbation  on  the  structure  relies  on  the  physical  constraint 
that  g(  r)  =0  for  r<d  when  the  pair  Interactions  contain  a  hard  core 
repulsion.  Once  the  thermodynamic  and  structural  properties  of  the 
trial  system  are  obtained  (In  some  degree  of  approximation),  It  Is 
then  quite  easy  to  relate  these  results  to  the  corresponding  proper¬ 
ties  of  the  fluid  In  which  the  pair  potential  Is  w( r)  rather  than 
wT(r).  This  Is  the  program  we  adopt  In  the  present  work. 

A  second  reason  for  separating  w(  r)  Into  a  harshly  repulsive 
reference  potential  plus  c  slowly  varying  perturbation  has  to  do  with 
the  convergence  of  perturbation  theories.  While  the  effects  of  a  quickly 
varying  perturbation  on  the  structure  of  a  fluid  can  be  very  large,6 
the  effects  of  a  slowly  varying  perturbation  are  small  by  comparison.' 

As  a  result,  when  It  Is  possible  to  separate  w( r)  In  such  a  way  that 
u(r)  Is  slowly  varying,  a  perturbation  theory  for  describing  the  effects 
produced  by  u( r)  can  be  rapidly  convergent.  It  Is  for  this  reason  that 
the  potential  separation  exhibited  in  Fig.  1  leads  to  a  first  order  per¬ 
turbation  theory4  that  Is  significantly  more  accurate  than  that  of  Barker 
and  Henderson6  who  used  the  negative  part  of  the  Lennard-Jones  potential 
(rather  than  the  attractive  part)  for  their  perturbation. 

The  notation  we  use  to  distinguish  the  properties  of  various  systems 
Is  the  following:  Properties  of  the  hard  sphere  fluid  are  denoted  by 
a  subscript  "d"  which  stands  for  the  diameter  of  the  spheres.  The 
properties  of  the  trial  system  are  identified  with  a  subscript  "T". 

Those  of  the  reference  system  are  denoted  with  a  "0".  Properties  of 
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the  fluid  of  interest  have  no  subscripts.  Thus,  for  example,  y(  k) 
denotes  the  structure  factor  for  the  system  with  pair  potential  w(  r) , 
while  Xtj(k)*  XT(  k)  and  Xo(  k)  denote  the  structure  factors  for  the 
hard  sphere,  trial  and  refere-ice  systems,  respectively. 

In  the  next  part  of  this  'ection,  we  establish  the  fundamental 
relationship  between  the  trial  system  and  the  actual  system  under 
consideration.  In  Part  C,  the  high  temperature  approximation  ( HTA) 
will  be  applied  to  the  trial  system,  and  the  relationship  wi 1 1  be 
used  to  obtain  the  corresponding  approximation  for  the  fluid  of  in¬ 
terest.  In  Part  D,  the  optimized  random  phase  approximation  (ORPA) 
will  be  used  in  the  same  way  to  obtain  more  accurate  results.  Some 
limitations  of  HTA  and  ORPA  are  discussed  in  Part  E. 
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B.  Relationship  between  Trial  System  and  Actual  System 

The  properties  of  a  system  with  realistic  repulsive  forces  can 
be  related  to  a  system  with  hard  core  repulsions  by  a  cluster  expan¬ 
sion.  The  expansion  developed  below  is  a  generalization  of  one  re¬ 
cently  presented  by  the  authors’. 

Note  that  Ct  is  a  functional  of  the  Mayer  function  f ( r)  =  exp(-pw)-l. 
The  corresponding  quantity  for  the  trial  system,  depends  In  the 

same  way  on  fT(r;d)  =  exp(-pw^.)-1.  These  cluster  functions  are  drawn 
In  Fig.  2a.  The  difference  between  them,  Af  =  f-fjr  Is  depicted  In 
Fig.  2b.  If  the  trial  system  diameter  d  Is  chosen  well,  the  distance 
over  which  Af(r)  Is  non-zero,  §d,  Is  small.  This  Immediately  suggests 
that  ^ and  can  be  related  by  the  following  functional  Taylor  series: 

&  =  ^T  +  I  r)  +  •  •  •  • 

The  subsequent  terms  are  higher  order  functionals  of  Af( r) .  The  effec¬ 
tive  expansion  parameter  In  this  series  Is  §,  which  can  be  regarded  to 
be  approximately  equal  to 

d  [  |Af(  r)  |  dr  . 

Jo 

By  performing  the  Indicated  functional  differentiation,  we  obtain 

d  =  dj  +  I  P*  J  yT(r)  Af(r)  dr  +  ...  ,  (9) 

where  yj(r)  Is  related  to  the  radial  distribution  function  In  the 
trial  system  b^r 


yT(r)  =  exp[+pwT(r)]  gT(r)  . 


(10) 
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From  Eq.  (9)  Jt  is  apparent  that  a  good  choice  of  d  5 s  that  value  for 
which  the  first  correction  to  ~  vanishes.  Thus,  d  is  chosen  to 
satisfy 

J  yT(r;d)  Af(r;d)  d^  =  0  .  (11) 

It  can  be  shown9  that  when  d  is  determined  by  Eq.  (11), 

*  0(5*)  .  (12) 

That  is,  Eq.  (11)  not  only  makes  the  term  of  order  5  in  Eq.  (9)  equal 
to  zero,  it  also  forces  further  corrections • to  be  of  order  5*  (and  not 
52).  From  the  usual  virlal  expansion  for  Cl*  it  Is  also  evident  that 
~  d)  is  exact  In  the  limit  of  small  p  when  d  is  chosen  by  Eq.  (11). 

A  similar  functional  Taylor  series  can  be  generated  to  relate 
h(r)  to  the  trial  system  correlation  function.  In  particular,  when  d 
is  chosen  by  Eq.  (11),  It  Is  found9  that 

y(r)  s  exp[+pw(r)]  g(r)  =  yT(r;d)  [1  +  0(§2)]  .  (13) 

As  a  result,  the  structure  factor  is  given  by 

X(k)  ■  XT(k;d)  +  p  J  dr  yT(r;d)  Af(  r)  exp(-ik-£)  +  0(§2)  (14) 

The  second  term  makes  appreciable  contributions  (in  the  form  of  a 
damping)  at  intermediate  and  large  wave-vectors.10 

Equations  (11) — ( 14)  provide  the  desired  relationship  between  the 
trial  system  and  the  system  in  which  the  pair  potential  is  w( r) .  We 
now  consider  two  approximate  solutions  to  the  problem  of  calculating 
properties  of  the  tr?al  system. 
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C.  High  Temperature  Approximation  for  the  Trial  System 

The  simplest  theory  for  the  trial  system  Is  obtained  by  neglect¬ 
ing  entirely  the  effects  of  the  perturbation  on  the  structure.  This 
approximation,  which  becomes  exact  in  the  limit  of  small  3,  is  called 
the  high  temperature  approximation  (HTA).  For  the  trial  system,  it 
is 

[  9^  r)l  HTA  =  9d^  0®) 

and 

^T^  HTA  =  ^d  “  2  Pp2  J  9d(r)  d~  *  (  l6) 

With  the  aid  of  Eqs.  (11) -(13)  we  obtain  from  (15)  and  (16)  the 
HTA  for  the  actual  system  with  the  pair  potential  w( r)  : 

9HTA(r)  =  exPl-0uo(r)]  YdJr)  =  9o(  r)  [  1  +  0(?1)]  (17) 

and 

^KTA  '  4.  '  2  J  Wr>  “<r>  dI  •  <l8> 

The  quantity  d0  is  the  solution  to  Eq,  (11)  when  the  HTA  is  used  for 
yT(  r;d)  ;  that  is 

J  Vdo(  r)  HU  dr  =  0  (19) 

where 

Afo  =  exp(-0uo)  -  exp(-pu.  )  . 

a  o 

Equations  ( 17) -( 19)  are  the  principal  equations  employed  by  the 
authors  in  their  recent  discussions  on  the  perturbation  theory  and 
structure  of  simple  liquids.1*4  At  high  densities  they  afford  an  ex¬ 
ceptionally  accurate  (as  well  as  simple)  theory  of  Leonard- Jones  liquids. 
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At  low  densities,  however,  the  structure  is  no  longer  dominated  entire 
ly  by  the  repulsive  forces,  and  the  attractions  must  be  accounted  for. 
Further,  it  is  probable  that  important  classes  of  systems  exist  (e.g. , 
those  in  which  the  attractive  forces  vary  rapidly  in  space)  for  which 
the  attractions  play  an  important  role  in  forming  the  liquid  structure 
even  at  high  densities.  Thus,  we  now  turn  to  the  problem  of  describ¬ 
ing  the  effects  of  the  perturbation  on  the  structure. 
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Dj _ Optimized  Random  Phase  Approximation  for  the  Trial  System 


Discussions  of  the  random  phase  approximation  (RPA)  and  the  op¬ 
timized  random  phase  approximation  (ORPA)  for  classical  fluids  have 
already  been  publ i shed. 5* 11 ~  1  3  Thus,  we  present  below  only  a  brief 
outline  of  the  basic  ideas,  together  with  the  formulas  necessary  for 
the  present  work. 

The  RPA  for  the  trial  system  describes  the  effects  of  the  per¬ 
turbation  on  the  structure  of  that  system  by  the  equation 

Xd(k) 

[xT(k)]RPA  -  i +p(k)Xd(k)  ’  ‘  (20) 

where 

p(k)  =  0p  ^(k)  =  p  J  djp  puT(r)  exp(-lk-£)  (21) 

This  approximation  can  be  expressed  in  terms  of  the  direct  correlation 
function,  c(r): 


[cT(r),RM  "  Cd(r>  *  *l“T(r>  ' 


(22) 


Thus,  it  is  the  first  order  approximation  In  the  perturbation  expan¬ 
sion  for  cT(r)  about  Cj(  r)  (where  the  ordering  parameter  is  the  strength 
of  the  perturbation  potential).  In  terms  of  the  radial  distribution 


function,  the  approximation  is 


[9T(r)]RPA  =  9d(r) 


sum  of  all  singly  connected  chains  In 
the  cluster  expansion  for  g^( r)  involv¬ 
ing  one  or  more  -pu^.(r)  bonds,  and  zero 
or  more  hd(  r)  bonds. 


=  9d(  r)  -  (  J  dk  */(k)  0v(k)  exp(  lk-_r)  , 


(25) 
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whe  re 

v(k)  =  uT( k)/[  1  +  p(k)  xd(k)]  .  (24) 

The  cluster  diagrams  described  in  Eq.  (23)  are  depicted  in  Fig.  3a. 

The  RPA  for  the  Helmholtz  free  energy  is 

-  2  Sp'  J  dZ  9d(r>  uT(r> 

+  \  (2tt)-3  J  dk  {p(k)xd(k)  -  £n[  1  +  p(  k)xd(  k)])  .  (25) 

Equations  (20),  (22)  or  (23)  are  obtained  by  taking  the  appropriate  func¬ 
tional  derivatives  of  (25)  [see  Eq.  (3)].  The  first  two  terms  on  the 
right  hand  side  of  Eq.  (25)  form  the  HTA  for  Q,-.  The  last  term  is 
the  sum  over  all  ring  graphs  In  the  cluster  expansion  for  in¬ 
volving  hj(  r)  bonds  and  two  or  more  -3u^.(r)  bonds.  These  diagrams  are 
depicted  in  Fig.  3b.  It  is  seen  that  the  RPA  can  also  be  appropriately 
called  the  "ring  picture  approximation". 

Since  the  RPA  represents  the  sum  of  a  class  of  graphs  in  the 
cluster  series,  systematic  corrections  to  the  approximation  can  be  ob¬ 
tained  b/  accounting  for  subsequent  terms  in  the  cluster  expansion. 

At  least  two  such  theories  have  been  presented:  the  yordered  expansion13 
and  the  mode  expansion. 5’ 1 1  ’ 1 2  The  former  utilizes  diagrammatic  methods 
explicitly,  while  the  latter  uses  a  collective  coordinate  formalism. 

Both  theories  give  corrections  to  the  RPA  in  the  form  of  an  infinite 
series  expressed  entirely  In  terms  of  reference  (i.e.  hard  sphere) 
fluid  correlation  functions  and  a  screened  potential  v( r) ,  which  Is  the 
Inverse  Fourier  transform  of  v(k)  defined  in  Eq.  (24). 

The  ORPA  is  a  method  for  enhancing  the  convergence  of  such  In¬ 
finite  series.  To  introduce  the  ORPA  we  note  that  the  short  ranged 
repulsive  interaction  In  the  trial  system  is  the  hard  core  potential 
Uj( r) .  Thus  Intermolecular  separations  smaller  than  d  are  physically 


-14- 


Inaccesslble.  As  a  result,  the  exact  evaluation  of  g^.(  r)  would  yield 
g^.{  r)  =  0  for  r<d.  However,  the  RPA  does  not  necessarily  obey  this 
condition  since  the  second  term  on  the  right  side  of  Eq.  (23)  may  be 
non- zero  for  r<d.  This  Improper  description  of  excluded  volume  effects 
is  related  to  the  fact  that  the  predictions  of  the  RPA  depend  on  the 
perturbation  in  this  physically  inaccessible  region12.  By  changing  u^.(  r) 
for  r<d  we  can  change  u^.(k)  (and  v(k))  for  any  value  of  k.  Such 
changes  obviously  alter  the  value  and  hence  the  accuracy  of 
and  the  form  of  [g^.(  r)]  RpA<  In  recent  work  on  model  electrolytes5*12 
it  was  found  that  by  choosing  u^.(  r)  for  r<d  such  that  [  g^.(  r)]  RpA  =  0 
for  r<d,  the  RPA  can  be  made  very  accurate.  Thus  we  are  motivated 
to  extend  the  ideas  used  In  the  electrolyte  calculations  to  the  prob¬ 
lem  under  consideration  here.  In  particular,  we  shall  pick  u(  r)  for 
r<d  so  that 

<Mr»RPA  -  -  (^1  -  0  •  <26> 

From  Eq.  (25)  It  is  seen  that  the  addition  of  this  physical  constraint 
to  the  RPA  Is  equivalent  to  demanding  that 

[6  aRpA/6u(r)]  =0  ,  r<d  ,  (27) 

where  aRpA  Is  the  random  phase  Integral 

aRPA  =  I  ( ?TT)"3  J  dJ5  {p(  k>Xd(  k)  -  tn[]  +  p(  k)Xd(  k)]}  ( 28) 

which  corrects  the  HTA  to  &j  in  the  RPA. 

In  summary,  the  ORPA  is  any  one  of  Eqs.  (20),  (22),  (23)  or  (25) 
supplemented  with  Eq.  (27).  Further  discussion  of  this  approximation 
is  contained  in  Reference  5. 
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E.  Limitations  of  HTA  and  ORPA 

Published  calculations  show  that  at  high  densities  the  HTA  gives 
a  very  accurate  theory  of  thermodynamic  properties  if  the  appropriate 
potential  separation  is  used.1 » 2 » 4  Further,  the  calculations  presented 
in  Section  IV  of  this  paper  show  that  the  ORPA  is  very  accurate  at  high 
and  moderate  densities  for  both  thermodynamic  and  structural  proper¬ 
ties.  This  state  of  affairs  does  not,  however,  persist  in  the  region 
of  low  densities.  In  particular,  the  HTA  and  ORPA  do  not  give  the 
exact  second  vtrtal  coefficient.  The  correct  second  virial  coefficient 
for  $  is 

B2  =  \  J  dr  [exp(-gu0  -  gu)  -  1]  . 

The  HTA  gives 

[Bz1HTA  =  2  I  I0-N  exp(-puo)  -  1]  , 

and  the  ORPA  gives 

tB^0RPA  =  2  I  d~  ^  1  '  pU  +  2  pZ  exP('Pu<>)  '  H  • 

The  error  in  the  ORPA  result  is  small  at  high  temperatures  but  can  be 
large  at  low  temperatures.  (The  corresponding  errors  in  the  low  density 
g( r)  are  similar.)  Away  from  the  critical  region,  the  errors  in  ORFA 
at  low  density  make  relatively  little  difference  in  the  measurable 
thermodynamic  properties,  since  for  temperatures  sufficiently  below  the 
critical  temperature  the  maximum  density  of  the  gas  is  very  small. 

The  HTA  and  ORPA  are  not  acceptable  theories  of  critical  phenomena. 
Neither  gives  a  free  energy  which  satisfies  thermodynamic  stability 
criteria.  Even  if  one  supplements  the  theoretical  free  energy  with  a 
Maxwell  construction,  a  purely  classical  coexistence  curve  is  obtained 
from  both  theories.14  Also  the  ORPA  (like  the  HTA)  is  inconsistent 


-16- 


In  the  sense  that  at  the  critical  point  (as  determined  by  the  free 
energy)  the  compressibility  as  determined  by  the  correlation  function 
Is  finite. 

One  pleasant  feature  of  ORPA  as  opposed  to  the  RPA  Is  that  the 
ORPA  does  not  contain  the  catastrophic  singularities  which  are  found 
In  the  RPA  and  frequently  (and  Incorrectly)  associated  with  phase 
transl  tlons. 1  5 
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III.  Calculational  Procedure  and  Summarizing  Equations 

The  first  step  when  applying  the  ORPA  is  the  determination  of  the 
optimized  form  of  u^.(  r)  for  r<d.  This  optimized  potential  is  the 
solution  to  the  variational  problem  stated  by  Eq.  (27).  We  find  this 
solution  numerically.  The  perturbation  in  the  unphysical  region, 
r<d,  is  expanded  in  a  series  of  basis  functions.  In  particular, 
u^.(  r)  Is  expressed  as 

-puT(r)  =  -pu(  r)  +  q  H(d-r)  -  (s/d)  G(  r-d) 


III 

+  H(d-r)  (1-r/d)2£  cn  fn#  (r/d) 


where  q,  s,  and  c^  are  constants,  H(x)  is  the  step  function 

H(x)  =  1,  x  >  0 
=  0,  x<0  , 

G(x)  is  the  function  whose  derivative  is  the  step  function 

G(x)  =  x,  x  >0 
=  0,  x<0  , 

and  (x)  =  Pn(2x-1),  where  Pp( y)  is  a  Legendre  polynomial.  The  Fourier 
transform  of  uT(  r)  thus  depends  on  the  values  of  q,  s,  c„,  . . .  ,  As  a 
result,  aRp^  also  depends  on  these  coefficients.  The  optimized  uT(  r) 
is  obtained  by  finding  the  set  of  coefficients  for  which  a^p^  is  station¬ 
ary  with  respect  to  small  changes  in  the  coefficients.  [We  note  that 
a^p^  is  a  positive  definite  functional  of  u^(  k) .  One  solution  of  Eq. 

(27)  corresponds  to  an  absolute  minimum,  and  it  is  apparently  this 
solution  that  we  find.  We  do  not  know  if  other  solutions  or  (27) 
exist.] 
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Wl th  the  coefficients  q,  s,  c0,  C|,  ...  determined  by  the  variational 
procedure,  the  RPA,  Eqs.  (20) -(25),  is  applied  directly  to  yield  the 
equilibrium  properties  of  the  trial  system.  In  particular,  Eq.  (23)  is 
applied  to  obtain  q^( r) .  In  order  to  relate  the  trial  system  proper¬ 
ties  to  the  actual  systems  properties,  we  also  need  y_  =  exp(+3wT)o  . 

For  r>d,  there  is  no  problem.  However,  the  ORPA  does  not  directly 
yield  Information  about  y^(  r)  for  r<d.  Fortunately,  we  need  y  (  r) 
for  values  of  r  close  tod  only.  Thus,  we  simply  extrapolate  into  the 
hard  core  using  the  fact  that  y^( r)  and  its  first  derivative  are  con¬ 
tinuous  functions  of  r.  The  y^  function  and  its  derivative  can  be 
evaluated  at  r-d  accurately  since  g^.(r)  -  c^(  r)  and  its  first  three 
derivatives  are  continuous  at  that  point.  Thus,  for  example, 


9j(d+)  =  cd(d+)  -  <=d(d-)  +  PuT(d”)  . 


The  direct  correlation  function  for  hard  spheres,  c.( r) ,  Is  known 
analytically,2  and  pu^d-)  can  be  evaluated  from  Eq.  (29). 

With  the  yT(r;d)  calculated  for  a  range  of  d  values,  we  then  cal¬ 
culate  the  d  which  satisfies  Eq.  (11).  It  Is  found  that  when  the  per¬ 
turbation  u(r)  Is  very  much  more  slowly  varying  than  u<>(  r) ,  the  value 
of  d  obtained  from  Eq.  (11)  is  essentially  the  same  as  d0  which  is  the 
solution  to  Eq.  (19).  With  the  value  of  d  determined,  Eqs.  (12),  (13) 
and  (14)  are  used  to  produce  the  final  ORPA  results  for  the  system  In 
which  the  pair  potential  is  w(r).  In  particular,  the  free  energy  is 


a. 


HTA 


+  Aa 


HTA 


+  a 


RPA 


(30) 


wh ere  Is  defined  by  Eqs.  (17)  and  (18)  (it  Is  a  function  of  the 
diameter  d0),  aRpA  is  the  random  phase  integral,  Eq.  (28),  evaluated 
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wfth  hard  sphere  diameter  d  and  the  opt imized  u^.(  r) ,  and 

4aHTA  =  ‘  2  6p!  J  dZ  <  VW  U(f)  ’  (3,) 

In  deriving  Eq.  (31),  we  make  use  of  Eq.  (19).  Since  for  many  systems, 
d  Is  very  close  In  value  to  do,  the  Is  frequently  negligible,  and 

aRpft(d)  can  be  replaced  by  aRRA(d0). 
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IV.  Application  to  Lennard-Jones  Fluid 

We  have  applied  the  procedure  described  in  the  previous  section 
to  the  Lennard-Jones  fluid,  which  is  a  hypothetical  fluid  of  spherical 
particles  which  interact  with  each  other  through  a  distance  dependent 
potential  of  the  following  form: 

w( r)  =  4e  [(o/r)12  -  (o/r)6]  .  (32) 

The  total  potential  energy  of  the  fluid  is  assumed  to  be  equal  to  the 
sum  of  all  the  pair  interactions. 

This  simple  model  is  known  to  be  inadequate  as  an  exact  descrip¬ 
tion  of  even  the  noble  gases,  because  of  the  restricted  form  of  the 
potential  and  the  neglect  of  three-body  forces.16  For  our  purposes 
its  major  virtue  is  that  Monte  Carlo  and  molecular  dynamics  calculations 
have  been  performed  for  this  model,  giving  information  about  the  thermo¬ 
dynamic  and  structural  properties.  Thus  for  checking  the  accuracy  of 
the  ORPA,  the  Lennard-Jones  fluid  Is  more  ideal  than  a  real  fluid  such 
as  argon,  for  which  we  do  not  know  the  intermolecular  interactions  so 
accurately. 

The  potential  separation  employed  in  our  calculations  is  the  one 
depicted  in  Fig.  1.  That  is, 


u0(r) 

«* 

u 

+ 

J 

n 

r  <  ro 

=  o 

r  >  ro 

(33a) 

u(r) 

=  "« 

r  <  ro 

=  w(  r) 

r>  r0  , 

(33b) 
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where  r0  =  2'^6  a  Is  the  position  of  the  minimum  in  w(  r) .  This  separation 
of  w(  r)  satisfies  all  the  conditions  desirable  for  reference  anJ  pertur¬ 
bation  interactions  (see  Section  II  A).  In  particular,  u0( r)  is  a  harsh 
repulsion,  and  u(  r)  is  slowly  varying. 

The  Lennard-Jones  parameters  c  and  o  can  be  used  as  the  units  of 
energy  and  length,  respectively;  all  the  thermodynamic  properties  can 
be  expressed  in  terms  of  the  dimensionless  density,  p  =  po3,  and  the 

dimensionless  temperature,  T*=(3e)~'.  In  these  units  the  critical  point 

*  #  # 

occurs  for  T£  =1.26  and  pc  =  .30,  and  the  triple  .joint  is  at  T(  =.70 

a 

and  p^  =.84.  The  accuracy  of  the  critical  point  parameters  is  uncer¬ 
tain.  1  7 

The  calculations  were  performed  for  a  wide  variety  of  states  for 
reduced  densities  f -om  about  .  1  to  1  and  reduced  temperatures  from  about 
.6  to  2.74.  For  each  state  the  hard  sphere  diameter,  do,  associated 
with  the  repulsive  part  of  the  potential  was  evaluated  [Eq.  (19)], 
followed  by  the  HTA  free  energy,  structure  factor,  and  pair  correlation 
function.  For  the  optimization  procedure  on  the  trial  system  we  used 
the  series  in  Eq.  (29)  truncated  after  the  P|  term.  [This  trial  function 
with  only  four  basis  functions  was  found  to  be  accurate  enough.  Ex- 
tending  the  series  up  to  ths  P5  term  had  no  appreciable  effect.] 

The  trial  system  hard  sphere  diameter,  d,  associated  with  the  repul¬ 
sive  part  of  the  potential  in  the  presence  of  attractions  was  evaluated 
according  to  Eq.  (11).  For  each  state  the  difference  between  d0  and  d 
was  found  to  be  negligible  (approximately  1-2x10“*  a).  The  optimized 
RPA  integral  provides  the  ORPA  contribution  to  the  free  energy,  and  the 


optimized  coefficients  were  used  to  calculate  the  ORPA  corrections  to 
the  structure  factor  and  pair  correlation  functions. 16 

These  calculations  require  a  knowledge  of  the  excess  free  energy 
and  the  pair  correlation  function  9j(  r)  f°r  t^ie  ^ar<^  sphere  fluid  for 


-22- 


range  of  densities.  We  used  the  results  of  Verlet  and  Wets2,  who  per¬ 
formed  Monte  Carlo  calculations  for  hard  spheres  and  presented  very 
convenient  analytic  forms  for  the  functions. 

The  thermodynamic  results  are  presented  in  Tables  I— III.  The 
free  energy  was  obtained  from  Eqs,  (18)  and  (30),  and  the  pressure 
and  energy  were  obtained  by  numerical  differentiations  of  the  free 
energy  with  respect  to  density  and  temperature,  respectively: 

■  -p  • 

and 

4E/N  .  .  , 

where  AE  stands  for  the  excess  (with  respect  to  the  ideal  gas)  internal 
energy.  These  results  are  compared  with  the  molecular  dynamics  re'sults 
of  Verlet19  and  the  Monte  Carlo  calculations  of  Levesque  and  Verlet, 
and  Hansen. 20  A  convenient  summary  of  much  of  the  available  computer 
data  is  found  in  Ref.  2. 

Excellent  agreement  between  ORPA  and  computer  results  are  obtained 
for  most  of  the  states  investigated.  For  example,  Verlet  estimated 
that  for  the  high  density,  low  temperature  states,  the  error  in  the 
molecular  dynamics  results  for  gp/p  may  reach  .05.  Table  II  shows 
that  for  most  of  the  high  density  states  the  differences  between  the  ORPA 
and  molecular  dynamics  results  are  equal  to  or  less  than  this  amount. 

The  only  serious  discrepancies  between  ORPA  and  Monte  Carlo  occur  for 
the  low  temperature  lav  density  gas,  l.e.  the  vapor  well  below  the 
critical  temperature. 

The  HTA  and  ORPA  structure  factors  and  pair  correlation  functions 
are  compared  with  molecular  dynamics  results21  In  Figs.  4-7.  Figure 
4  shows  the  structure  factor  for  p  =.844,  T  =.723,  a  state  near 
the  triple  point.  The  major  discrepancy  between  the  HTA  and  molecular 
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dynamics  results  occur  only  in  the  height  of  the  first  peak.  The  ORPA 
raises  the  peak,  eliminating  much  of  the  discrepancy,  and  makes  other 
changes  which  are  not  visible  on  the  scale  of  this  graph.  Fig.  5 
shows  the  pair  correlation  function  for  a  nearby  state.  The  HTA  is 
remarkably  accurate  but  peaks  slightly  too  high  and  for  too  small  an 
intermolecular  distance.  The  ORPA  improves  the  agreement  with  computer 
experiments. 

These  two  figures  illustrate  the  accuracy  of  the  hypothesis  that 
at  high  density  repulsive  forces  dominate  the  structure  of  a  simple 
liquid.1  The  HTA  result  is  not  accurate  enough  to  give  good  pressures 
via  the  virial  equation. 2  This,  however,  is  an  extremely  sensitive 
test  of  a  correlation  function.  When  directly  compared  with  the  re¬ 
sults  of  the  computer  experiment,  the  HTA  can  be  seen  to  account 
quite  well  for  the  important  features  of  g( r) . 

Figures  6  and  7  show  the  corresponding  results  for  states  at 
about  critical  temperature  but  at  a  density  somewhat  higher  than  criti¬ 
cal.  For  g( r)  the  HTA  is  only  qualitatively  correct.  For  the  struc¬ 
ture  factor,  the  HTA  is  again  remarkably  g<5od  except  for  small  k.  The 
ORPA  to  a  large  extent  corrects  these  discrepancies.  For  the  structure 
factor,  the  ORPA  results  essentially  coincide  with  the  molecular  dynamics 
results,  except  at  small  k.  (The  small  k  errors  in  the  ORPA  are  a  re¬ 
flection  of  the  fact  that  the  ORPA  does  not  properly  account  for  the 
long  range  correlations  in  the  critical  region.)  These  results  con¬ 
firm  the  conclusions  of  our  previous  work1  about  the  density  and  k 
dependence  of  the  effect  of  repulsive  forces  on  the  structure  of  liquids  22 

The  ORPA  results  also  agree  quite  closely  with  the  second  order 


perturbation  theory  of  Barker  and  Henderson.6  For  example,  along 
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the  isotherm  T#  =  1.15,  the  free  energies  (-^?/p)  of  the  two  theories 
differ  by  at  most  .02  and  the  compressibility  factors  j3p/p  by  at  most 
.06. 

A  major  difference  between  the  two  theories  is  that  at  high 
densities  the  free  energy  series  whose  first  two  terms  are  the  HTA 
and  ORPA,  converge  much  more  rapidly  than  the  BH  perturbation  theory. 

The  HTA  is  already  essentially  correct  at  high  densities,  and  the  ORPA 
makes  a  very  small  correction  (see  for  example  Table  III).  By  com¬ 
parison  trie  BH  first  order  theory  Is  not  nearly  so  accurate*  and  the 
BH  second  order  correction  Is  about  20  times,  larger  than  the  RPA  correc¬ 
tion.  A  similar  difference  in  the  convergence  of  g(  r)  and  x( k)  exists 
for  the  two  theories.  The  origin  of  these  differences  is  the  methods 
used  to  separate  the  potential  into  reference  and  perturbation  parts. 

Barker  and  Henderson  choose  the  positive  p?rt  of  the  Lennard-Jones  po¬ 
tential  as  their  reference  Interaction.  The  negative  part  is  taken  to 
be  the  perturbation.  However,  the  negative  part  of  the  Lennard-Jones 
potential  i s  not  slowly  varying,  and  thus  it  makes  large  contributions 
to  the  structure  at  all  densities.  The  separation  Into  purely  repul¬ 
sive  and  purely  attractive  parts  which  we  use  is  more  useful,  since  at 
high  densities  the  purely  attractive  part  (which  is  slowly  varying) 
has  only  a  very  small  effect  on  the  structure. 

A  second  major  difference  is  that  in  the  BH  theory,  it  Is  necessary 
to  have  Information  about  the  three-  and  four-particle  correlation  func¬ 
tions  in  the  hard  sphere  fluid  in  order  to  calculate  the  effect  of  the 
perturbation  potential  on  the  structure.  On  the  other  hand,  the  ORPA 
requires  only  the  more  easily  obtained  two-particle  g( r)  for  hard  spheres.25 
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V.  Discussion 

In  this  section,  we  summarize  the  principal  conclusions  of  this  work 
and  discuss  the  advantages  associated  with  using  the  ORPA  to  study 
classical  fluids. 

Principal  conclusions.  The  structure  of  dense  simple  classical 
fluids  is  determined  mostly  by  the  rapidly  varying  parts  of  the  inter- 
molecular  potential.  If  all  the  rapidly  varying  parts  are  short  ranged  and 
purely  repulsive,  a  quickly  convergent  theory  of  the  equilibrium  proper¬ 
ties  of  the  dense  fluid  can  be  constructed.  In  this  theory  the  effect  of 
the  harshly  varying  repulsions  is  related  to  the  properties  of  the  hard 
sphere  fluid,  and  the  effect  of  the  remaining  perturbation  on  the  structure 
is  calculated  from  the  optimized  random  phase  approximation.  The  resulting 
theory,  in  addition  to  being  rapidly  convergent  at  high  densities  is  ac¬ 
curate  and  useful  for  the  entire  single-phase  fluid  region  of  the  phase 
diagram,  excepting  the  critical  region  and  the  very  low  temperature  vapor. 

At  high  density  the  effect  of  the  perturbation  on  the  structure  can  be 
neglected  completely,  giving  a  very  accurate  thermodynamic  theory  and  a 
fairly  accurate  and  extremely  simple  .tpproximat  ion  for  the  pair  correla¬ 
tion  function,  g(  r)  ~  g  0(  r) . 

If  the  atti active  part  of  the  potential  has  spatial  variations  much 
greater  than  those  of  the  attractive  part  of  the  Lennard-Jones  potential 
(or  if  there  are  some  long  ranged  rapidly  varying  repulsive  forces  or  oscil¬ 
lations  in  the  potential  such  as  are  postulated  for  liquid  metals),  then 
these  parts  of  the  potential  can  influence  the  structure  significantly. 

In  particular,  according  to  the  ORPA  Eq.  (20)  for  X(  k) ,  if  the  perturbation 
has  appreciable  negative  Fourier  components  near  k  =  2t/(J,  where  Xo(  k)  is 
large,  then  the  perturbation  wi 1 1  affect  the  structure  at  high  density.24 
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Advantages  of  the  optimized  random  phase  approximation.  Computations 
involving  the  optimized  random  phase  approximat  *n  requires  detailed  but 
available2  information  about  the  free  energy  and  pair  correlation  function 
of  the  hard  sphere  fluid.  The  accuracy  obtained  with  just  this  informa¬ 
tion  is  comparable  to  that  obtained  when  one  carries  the  Barker-Henderson 
perturbation  series  for  Ci  t<>  second  order  in  the  perturbation.6  However, 
the  second  order  term  in  the  perturbation  series  requires  knowledge  about 
of  three-  and  four-particle  correlations  as  well  as  the  two-particle  cor¬ 
relation  in  the  hard  sphere  fluid.  With  the  same  information  needed  to 
carry  the  straightforward  perturbation  series  to  second  order,  one  can 
use  either  the  mode  expansion5’ 1 '» 12  or  y-ordered  expansion13  to  calculate 
correct  Ions  to  the  ORPA.  For  some  fluids  the  corrections  may  be  necessary 
to  obtain  quantitative  results  when  the  perturbation  Is  significantly 
different  from  the  Lennard-Jones  perturbation. 

Another  advantage  of  the  ORPA  Is  that  It,  unlike  the  weak  interaction 
perturbation  expansion,25  is  particularly  well  suited  for  describing  long 
ranged  Interactions,  such  as  the  Coulomb  potent ial . 5*  1 '»  12  Finally,  the 
ring  and  chain  summations  which  give  the  random  phase  approximation  can 
also  be  performed  for  a  certain  class  of  highly  anisotropic  molecular  po¬ 
tentials.7’26  The  generalized  RPA  formula  for  these  angular  dependent 
forces  are  only  slightly  more  complicated  than  those  for  spherically  sym¬ 
metric  forces.  Thus,  techniques  similar  to  those  discussed  in  this  paper 
can  be  applied  to  many  types  of  fluids  which  are  more  Interesting  and  com¬ 
plex  than  single  liquids. 
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Jable  I.  Excess  free  energy  for  Leonard- lonf»c  f  i„;j  r  1  _ 
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— ?ble  U.  Pressure  and 
several  isotherms. 
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Table  II.  Pressure  and  internal  energy  for  Lennard-Jones  fluid  along 
several  isotherms. 
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4.21 

9.58 

3.75 

-AE/Ne 


0RPAb 

MCC 

.  85 

1.15 

5.06 

5.07 

5.77 

5.78 

6.03 

6.04 

.74 

.86 

1.47 

1.55 

2.  14 

2.24 

2.77 

2.  85 

3.43 

3.47 

4.11 

4.  14 

4.46 

4.45 

5.  12 

5.  13 

5.68 

5.67 

5.97 

5.96 

.71 

.78 

1.39 

1.51 

2.04 

2.09 

2.69 

2.75 

3.35 

3.37 

3.68 

3.70 

4.67 

4.68 

5.26 

5.25 

5.68 

5.66 

5.81 

5.71 

.60 

.61 

1.19 

1.21 

1.77 

1.78 

2.36 

2.37 

3.21 

3.21 

3.94 

3.90 

4.28 

4.28 

4.40 

4.41 

4.22 

4.18 

3.76 

3.  80 

a. 

b. 

c. 


Obtained  from  numerical  differentiation 
Obtained  from  numerical  differentiation 
Monte  Carlo  computer  experiment  results 


°f,SHTA- 

of  <W 

(see  Ref. 


2). 


d.  For  this  thermodynamic  state,  pd!>.93.  The  analytic  gd(  r) 
and  #d  used  in  the  theoretical  calculations  presented  here 
are  reliable  only  for  pd!<.93. 


33- 


Table  III.  Pressure  and  internal  energy  for  Lennard-Jones  fluid  along 
several  high  . ensity  Isochores. 


Pp/P 


p* 

T* 

HTAa 

0RPAb 

.88 

1.095 

3.44 

3.47 

.88 

.940 

2.  82 

2.  85 

.88 

.591 

-.32 

-.26 

:  85 

2.889 

4.38 

4.39 

.85 

2.202 

4.22 

4.23 

.85 

1.214 

3.04 

3.07 

.85 

1.128 

2.78 

2.  81 

.85 

.88 

1.64 

1.69 

.85 

.782 

.93 

.98 

.85 

.76 

.74 

.79 

.85 

.658 

-.35 

-.29 

.85 

.591 

-1.32 

-1.24 

.75 

2.849 

3.  10 

3.11 

.75 

1.304 

1.55 

1.59 

.75 

1.069 

.75 

.81 

.75 

1.071 

.76 

.82 

.75 

.881 

-.27 

-.  18 

.75 

.827 

-.67 

-.57 

.65 

2.557 

2.11 

2.12 

.65 

1.585 

1.19 

1.24 

.65 

J.036 

-.  32 

-.21 

.65 

.90 

-1.04 

-.  89 

-AE/Ne 


MDC 

HTAa 

0RPAb 

MDC 

3.48 

5.86 

5.  87 

5.85 

2.77 

6.03 

6.04 

6.04 

-.18 

6.44 

6.46 

6.53 

4.36 

4.25 

4.26 

4.25 

4.20 

4.77 

4.78 

4.76 

3.06 

5.60 

5.62 

5.60 

2.78 

5.69 

5.70 

5.69 

1.64 

5.93 

5.95 

5.94 

.98 

6.04 

6.05 

6.04 

.82 

6.06 

6.08 

6.  07 

-.20 

6.  17 

6.19 

6.39 

-1.20 

6.25 

6.27 

6.46 

3.10 

4.05 

4.07 

4.07 

1.61 

4.99 

5.01 

5.02 

.90 

5.  15 

5.  18 

5.19 

.89 

5.15 

5.  17 

5.  17 

-.12 

5.29 

5.32 

5. 31 

-.54 

5.33 

5.36 

5.38 

2.14 

3.77 

3.80 

3.78 

1.25 

4.20 

4.24 

4.23 

-.11 

4.47 

4.52 

4.52 

-.74 

4.54 

4.60 

4.61 

a.  Obtained  from  numerical  differentiation  of 

b.  Obtained  from  numerical  differentiation  of  £?0Rp^- 

c.  Molecular  Dynamics  computer  experiment  results  (see  Ref.  2). 
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Figure  Captions 


Fig.  1.  Intermol ecular  interactions.  The  Lennard-Jones  potential  is 

w(  r)  =  Uo(r)  +  u(  r) .  The  trial  system  interaction  is  the  hard 
core  potential,  u r) ,  plus  u(  r) . 

Fig.  2.  Mayer  cluster  functions,  a.  f  and  are  the  cluster  functions 
for  the  actual  system  and  trial  system,  respectively,  b.  The 
difference  Af  =  f-fy.  The  distance  ro  is  the  range  of  the  re¬ 
pulsive  potential  uo(r),  while  d  is  the  hard  sphere  diameter  in 
the  tr i al  system. 

Fig.  3.  Cluster  diagrams  in  the  random  phase  approximation,  a.  The 

radial  distribution  function,  b.  The  free  energy.  Solid  lines 
denote  -j3uy(  r)  bonds,  dashed  lines  are  h^(  r)  =  g^rj-l  bonds, 
root  points  are  open  circles,  and  density  field  points  are  filled 

T 

(i.e. ,  the  first  two  terms  on  the  right  hand  side  of  Eq.  (25)). 

Fig.  4.  High  density  structure  factor  for  the  Lennard-Jones  fluid.  The 
circles  denote  the  molecular  dynamics  results  (Ref.  21);  the 
dashed  line  represents  the  high  temperature  approximation  (i.e,, 
X(  k)  -  Xo(k));  the  solid  line  is  the  optimized  random  phase 
approximat ion. 


circles. 


^T^  HTA 


is  the  high  temperature  approximation  for  & 


Fig.  5.  High  density  pair  correlation  function  for  the  Lennard-Jones 

fluid.  The  circles  denote  the  molecular  dynamics  results  (Ref. 
21);  the  dashed  line  represents  the  high  temperature  approximation 
(i.e.,  g(  r)  ~  g  o(  r)) ;  the  sol  id  line  is  the  optimized  random 
phase  approximation. 
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Fig.  6.  Moderate  density  structure  factor  for  the  Lennard-Jones  fluid. 

The  circles  denote  the  molecular  dynamics  results  (Ref.  21);  the 
dashed  line  represents  the  high  temperature  approximation  (i.e. , 
X(k)  ~  Xo(k));  the  solid  line  is  the  opt imized  random  phase 
approximation. 

Fig.  7.  Moderate  density  pair  correlation  function  for  the  Lennard-Jones 
fluid.  The  circles  denote  the  mol ecular  dynamics  results  (Ref. 
21);  the  dashed  line  represents  the  high  temperature  approximation 
(  i.e. ,  g(  r)  ~  g0( •'))  ;  the  sol  id  line  is  the  opt  imized  random 
phase  approximation. 
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